Streaking and Wigner time delays in photoemission from atoms and surfaces 
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Streaked photoemission metrology allows the observation of an apparent relative time delay be- 
tween the detection of photoelectrons from different initial electronic states. This relative delay 
is obtained by recording the photoelectron yield as a function of the delay between an ionizing 
ultrashort extended ultraviolet (XUV) pulse and a streaking infrared (IR) pulse. Theoretically, 
photoemission delays can be defined based on i) the phase shift the photoelectron wavefunction 
accumulates during the release and propagation of the photoelectron ("Wigner delay") and, alter- 
natively, ii) the streaking trace in the calculated photoemission spectrum ("streaking delay"). We 
investigate the relation between Wigner and streaking delays in the photoemission from atomic and 
solid-surface targets. For solid targets and assuming a vanishing IR-skin depth, both Wigner and 
streaking delays can be interpreted as an average propagation time needed by photoelectrons to 
reach the surface, while the two delays differ for non-vanishing skin depths. For atomic targets, the 
difference between Wigner and streaking delays depends on the range of the ionic potential. 

PACS numbers: 42.65.Re, 79.60.-i, 
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I. INTRODUCTION 

Streaked photoemission spectroscopy is increasingly 
applied to resolve ultra-fast electronic processes at the 
natural timescale (« 1 atomic unit = 2.4 x 10~ 17 s = 24 
attoseconds (as)) of the motion of valence electrons in 
matter. Streaking metrology uses ultra-short pulses of 
extreme ultraviolet (XUV) radiation to emit electrons 
into the electric field of a delayed infrared (IR) laser 
pulse. The XUV and IR pulses in this pump-probe setup 
are phase coherent, and the yield of emitted photoelec- 
trons is recorded as a function of the delay At between 
the two pulses [l|. The resulting energy-resolved pho- 
toemission spectra show stripes that oscillate in At with 
the period of the IR-laser electric field. These streak- 
ing traces occur in distinct photoelectron kinetic energy 
intervals that are determined by the spectral width of 
the XUV pulse, the IR-laser intensity, and the density of 
states of the target. For photoemission out of energeti- 
cally resolved discrete atomic levels, streaking traces can 
be related to a given initial state, while for photoemission 
from electronic states in solids they are modulated by the 
density of states within a given band of occupied states. 
The analysis of streaked photoemission spectra proceeds 
by fitting the center of energy (COE) of a given streak- 
ing trace as a function of At to a sine function with an 
adjustable phase, thereby mapping energy shifts induced 
by the IR laser onto a time delay between the apparent 
release of the photoelectron and its arrival of the XUV 
pulse @- 

Streaked photoelectron spectra from localized states 
of atomic targets Q and both, delocalized conduction 
band and localized core-level bands of solids Q have re- 
cently been recorded, leading to an ongoing debate about 
i) the interpretation of the deduced photoemission time 
delays and ii) the possibility of distinguishing delay con- 
tributions from the primary XUV-photorelease process 
and subsequent photoelectron propagation in the ionic 



potential of the target and the streaking IR-laser elec- 
tric field Applying streaking metrology to neon 
atoms, Schultzc et al. [3j measured a relative photoe- 
mission streaking delay of At$ — 21 ± 5 as for the 
release of electrons from 2p orbitals relative to emis- 
sion from 2s orbitals. The authors theoretically ana- 
lyzed their measured relative delay in terms of the rel- 
ative Wigner delay Atvk(c) that is given by the energy 
derivative of the spectral phase of the calculated pho- 
toelectron wave function |g, [l(J and averaged Arvy(e) 
over the spectral profile of the XUV pulse. The differ- 
ence between their calculated averaged relative Wigner 
delays Atw for emission from the 2s and 2p orbitals 
did not exceed Ar(2p — 2s) = rw(2p) — Tw(2s) = 6.4 
as, even for calculations that included electronic correla- 
tion in neon at the multiconfiguration-Hartree-Fock level. 
The authors linked the mismatch between their measured 
relative streaking and calculated Wigner delays to the 
extreme sensitivity of the photoelectron wavefunction's 
spectral phase to electronic correlation effects in multi- 
electron atoms. 

The same experiment [|| was subsequently analyzed by 
Khcifcts and Ivanov [f| based on numerical solutions of 
the time-dependent Schrodingcr equation (TDSE) for a 
single active electron moving in the Hartree-Fock poten- 
tial of a Ne + ion and, in a separate approach, by including 
electronic correlation effects to some extent by numeri- 
cally solving a set of coupled equations in a random-phase 
approximation-with-exchange model. These calculations 
reproduce only less than one half of the measured relative 
delay of 21 as, and the authors speculated that the much 
larger observed relative delay might not be due solely 
to the XUV- induced release process, even if electronic 
correlation effects were accurately accounted for. The 
measured relative streaking delay might thus include sig- 
nificant contributions from the photoelectron's interac- 
tion with the streaking IR-laser electric field. Indeed, the 
single-electron TSDE calculations by Ivanov [?J showed 
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that the IR electric field has a considerable influence on 
the Wigner delay for photocmission. The question then 
arises whether this IR-dresscd Wigner delay can be used 
to interpret the measured relative streaking delay. As 
shown in our previous investigation of time-resolved pho- 
toemission from a one-dimensional (ID) model hydrogen 
atom Q , the streaking delay is independent of the streak- 
ing laser field intensity. This was confirmed recently in 
a full-dimensionality calculation for atoms by Nagele et 
al. @. 

Investigating photocmission from a tungsten surface, 
Cavalieri et al. [J] have measured a relative streaking de- 
layof At 5 (OB-4/) = T S (CB)-T S {if) = 110±70asfor 
electrons emitted from 4f core levels relative to electrons 
released from the conduction band. This relative delay 
was interpreted as the delayed onset of IR streaking, i.e., 
as the difference in time needed by 4f and conduction 
band electrons to reach the surface [3, [n|. However, as 
we will argue in Sec. IIII Bl below, this interpretation is 
only valid under the assumption that the streaking IR 
field is fully screened inside the solid. We will also show 
in this work that Wigner and streaking delays become 
identical only in this limit of a sudden onset of IR streak- 
ing at the surface. We will show that the streaking delay 
sensitively depends on the IR-skin depth 8l- Therefore, 
the intuitive interpretation of relative streaking delays 
in terms of an effective photoelectron path length inside 
the solid becomes questionable for realistic values of 6l, 
depending on how exactly the IR electric field becomes 
screened in the solid. 

While only relative streaking delays can be deduced 
from measured photocmission spectra, Wigner time de- 
lays are conveniently derived from calculated photoelec- 
tron wave functions. Wigner and streaking delays in 
time-resolved atomic photocmission were examined re- 
cently 0, Q , and the nature of delays within the general 
context of scattering, decay, and photo- and particle- 
induced emission processes in atomic, nuclear, other 
branches of physics has been discussed by theorists for 
more than half a century [ij [To| (for a recent review 
see [l2j ■) In this work, we investigate the relation be- 
tween Wigner and streaking delays for photoemission 
from atoms in the gas phase and solid surfaces. In 
Sec. im we present the underlying theoretical models 
and our schemes for calculating time delays in photoe- 
mission. In Sec. IIIIl we compare and discuss our numer- 
ical results of time-resolved photoemission spectra from 
atoms (Sec. IIIIA ) and core-level and conduction bands 
of solid targets (Sec. IIIIB ). In particular, we investigate 
the dependence of the corresponding time delays on the 
XUV-photon energy, the effective range z c of the atomic 
model potential, initial state, IR-skin depth, and position 
of the Fermi level. Our conclusions follow in Sec. IV. 
Unless indicated otherwise, we use atomic units (a.u.) 
throughout this work. 



II. DEFINITION AND COMPUTATION OF 
WIGNER AND STREAKING TIME DELAYS 

The essence of the time delay introduced by Wigner 
and Smith 0, [13] can be understood for the elementary 
example of potential scattering in one spatial dimension, 
by representing the projectile as an incident wave packet 

61>in(z,t) = Jdkake 1 ^-^ (1) 

in terms of a superposition with amplitudes ak of plane 
waves with momenta k centered about k c and energies 
£fc. Scattering subject to the projectile-target-interaction 
potential V of finite range results in the outgoing wave 

SiP out (z,t) = J dka k e^e lkz -* ekt (2) 

for which each spectral component is phase shifted by tpk 
relative to the corresponding component of the incident 
wave. Depending on the values of the scattering phase 
shifts tfk , wave fronts and center of the scattered wave are 
shifted relative to the incident wave. The phase shifts (pk 
thus quantify the effect of V on Sipi n . Depending on the 
nature of V, wave fronts of all plane- wave components 
and the center and crest of the scattered wave packet 
may appear behind or ahead of the corresponding terms 
of the incident wave packet. The scattered wave packet 
can thus be characterized by a positive or negative delay 
time r, depending on whether its wave fronts or center 
are detected after or before they would be detected in 
the absence of V, respectively. More precisely, the phase 
shifts of individual traveling plane wave components lead 
to spectral delays 

Tw(£k) = — , (3) 

which, evaluated at the spectral center e c = k\l2 of the 
incident wave packet, define the Wigner delay [9. Hoj 

Tw=r w (e c ). (4) 

An alternative method for assessing the delay of the 
scattered relative to the incident wave packet is given in 
terms of the expectation values for the position, 

/>oo 

(«)(*)= / dzz\8^ out {z,t)\ 2 , (5) 
Jo 

and velocity, 

/>oo 

(v) = / dkkScj> out (k,t)\ 2 , (6) 
Jo 

of the scattered wave packet at a sufficiently large time 
t > T after the interaction according to [HI, [l]| 

<*> = <«>(*-*#)> (7) 
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FIG. 1: (Color online) Interpretation of the Wigner delay in 
the photoionization of atoms. The solid line shows schemati- 
cally the expectation values for the position (z) of the photo- 
electron wave packet. The photoemission Wigner delay Tyy is 
determined by a straight-line extrapolation of (z) according 
to Eq. (7) (dashed line). The XUV pulse is centered at t = 0. 



where 8<j) ou t(k,t) is the Fourier transformation of 
StjJout(z,t). The time T is is chosen so that V((z)(t > 
T)) ru and all spectral phases <p k and (v) of the outgo- 
ing wave packet remain time- independent to a very good 
approximation for t > T . This fitting procedure can be 
understood classically by identifying the expectation val- 
ues for position and velocity in Eq. ([7]) with the motion 
of a point particle. We refer to both tI v and t^, as 
"Wigner delay" since it can be shown |12| that they are 
closely related by the expression 



~w 



dk\6(j>out(k,t)\ 2 Tw(ek)- 



(8) 



Indeed, our numerical results in Sec. Mil below confirm 
that Tyy and t^J are almost identically. 

Both definitions of the Wigner time delay, and 
Tyy , have been used to characterize the scattering pro- 
cesses fl2l ] and can also be applied to determine delays in 
photoemission. They were recently used to interpret the 
relative delay in the IR streaked XUV spectra from the 
2s and 2p shells in neon [3, 0, 0] • We thus believe that it 
is important to carefully investigate the relation between 
the Wigner and streaking delays. 

We calculate the two Wigner delays for photoioniza- 
tion by solving the TDSE for photoelectron wave packets 
5tp(z,t) emitted from an initial state ipi(z,t) [f|, 



d 



' 2 



V(z) 



6ip(z,t) 



zE x (t)^i(z,t), 



(9) 



where p = id/dz is the momentum operator and 
il>i(z,t) = e^ 6B ^^pi(z) the stationary initial state with 
binding energy eb- For convenience we drop the sub- 
script "out" and designate the outgoing photoelectron 
wave simply as Stp(z,t). We consider photoemission by 
an attosecond XUV pulse from either the ground state of 



a ID model atom or from the energetically lowest bands 
of occupied initial states in the periodic ID model poten- 
tial of a solid. For the atomic case, we choose tpi(z,t) to 
be the ground state of the model atom with a potential 
V(z). For photoemission from a surface, ipi(z,t) desig- 
nates Bloch waves from within a given band with Block 
momenta k%. 

We represent the coupling of the XUV pulse electric 
field Ex to the electron in the dipole-length form and 
assume a Gaussian pulse profile, 



E x {t) 



-21n2(t/T X ) 2 sin( ^ 



(10) 



We assume a pulse duration tx = 300 as and a variable 
central frequency uix- 

We propagate the photoelectron wave packet from 
—4 fs to 8 fs on a spatial grid that extends over 32,000 a.u. 
and calculate the Wigner delays r™- according to ((4]) 
and ([7]) . Rewriting the photoelectron wave packet at the 
large time T > t x as [1, 



8Tp(z,t = T)~ dk^ k {z)d k Ex{e k -e B )e 



-ie k T 



(11) 



where ipk(z) is a continuum eigenstate in the potential 
V(z) with energy e k , d k = (^feM^i) the dipole matrix 
element, and Ex{uj) — J dtEx{t)e lult the spectrum of 
Ex{t), demonstrates that the Wigner delays depend on 
both, the electron release during the dipole coupling of 
the XUV electric field and the propagation of the photo- 
electron in the continuum. 

In order to calculate streaking delays from the streaked 
photoemission spectra, the influence of the streaking IR 
laser field on the active electron need to be investigated. 
The IR streaking effect on the release and propagation 
of the photoelectron is included by replacing p with 
p + A L (z,t - At) in Eq. (J9]). The influence of the IR- 
laser on the initial state can be included by numerically 
propagating the initial state in the IR-laser electric field 
according to 0, 

i^Piiz, t) = | i [p + A L (z, t - At)} 2 + V{z)\ iPi{ Z> t), 

(12) 

where At is the delay between the centers of the XUV 
and IR pulses, and the convention is used that At > 
corresponds to the XUV pulse preceding the IR pulse. 
We model the vector potential of the IR-laser pulse as 

A L (t) = A) sin 2 (nt/T L ) cos [ujl (t — tl/2)] (13) 

for < t < tl and set to otherwise. As pulse 
parameters, we choose the central photon energy Hlul = 
1.57 eV (corresponding to a wavelength of Al = 800 nm), 
peak intensity II = = 5x 10 n W/cm 2 , and pulse 

length tl — 8 fs. 

Assuming a free-electron dispersion (e = k 2 /2), the 
energy-differential photoemission probability is given by 



P(e, At) = i \6c/)(k,oo;At)\' 
k 



(14) 
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where 54>(k, oo; At) is the Fourier transform of 6ip(z, t — ¥ 
oo; At). The XUV-IR delay-dependent COE for a given 
streaking trace is 0] 



EcoE(At) 



1 



2P tot (At) 
with the total emission probability 

Ptot(At) = { dk\S<j>{k,oo;At)f 



dk \k 5(/)(k, oo; At)\ 2 , (15) 



(16) 



After calculating EcoE(At) for a range of XUV-IR de- 
lays — tx/2 < At < tl/2, we obtain the streaking delay 
ts relative to Ai bv fitting the parameters a, b, and ts 
to the expression [J, Q 



E C OE{At) =a + bA L (At - ts) 



(17) 



For XUV photocmission from solids, Eqs. (fT4"|) -(fT7 j) re- 
main valid. In this case, the initial states in Eq. ^ are in- 
dividual Bloch waves with momenta in the first Brillouin 
zone of either core-level or conduction band. We will 
show in Sec. IIII Bl below how to calculate band-averaged 
results. 



III. NUMERICAL RESULTS 
A. One-dimensional model hydrogen atom 

In this section, we discuss our numerical results for 
Wigner and streaking delays for XUV photoemission 
from the ground state with binding energy eb = 13.6 eV 
of the soft-core Coulomb potential 



V(z) = V c (z) = -l/Vz 2 + 2. 



(18) 



In the calculation of the Wigner delay r^, the direct nu- 
merical determination of the phase of the photoelectron 
wave packet according to tp k = \n5(f>(k,T)/\5(f>(k,T) is 
inaccurate or impossible for values of k where \S(j){k,T)\ 
is extremely small. Furthermore, e c is difficult to deter- 
mine from S4>(k, T) [see Fig. [2] (a)]. In order to overcome 
these two difficulties, we fit the real part Ke[5(f>(k, T)} of 
the calculated photoelectron wave packet to the function 



with 



/(e) = A e - 21n2 [( £ -^)/Ae] 2 cos (^ fc ) (ig) 



ifk = a(e-e c ) + /3(e - e c ) 2 + 7, (20) 



and determine (by least-squares fit over a large range 
of e = k 2 /2 values) the parameters e c , Ae, a = t-^,., 
(3 and 7. As shown in Fig [5] for photoionization with 
huix = 25 and 50 eV, the spectra of the calculated 
(by numerically solving the TDSE) and fitted photoelec- 
tron wave packets are in excellent agreement. Since ac- 
cording to Eq. ([TT]) the XUV pulse spectral profile is 
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FIG. 2: (Color online) Real part of the calculated and the 
fitted momentum-space photoelectron wave packets for ion- 
ization from ID model hydrogen atoms by XUV pulses with 
energies of (a) 25 eV and (b) 50 eV. 

imprinted on the photoelectron wave packet, we find 
that the fitted values Ae{hwx = 25 eV)=5.96 eV and 
Ae{hwx = 50 cV)=6.11 eV are close to the spectral 
width HAujx = 6.08 eV of the XUV pulse. 

Figure [3] shows the Wigner delays t^ 11 for the ID 
model hydrogen atom, calculated according to Eqs. ([7| 
and (fT9| - ((20l) . in comparison with the streaking delay ts, 
calculated as described in Sec. Ui. All delays are neg- 
ative. The two Wigner delays are almost identical and 
their absolute values are much larger than the streaking 
delay. Having established that the two Wigner delays al- 
most coincide for all parameters considered in this work, 
we only show results for r^J from now on and drop the 
superscript "II" for convenience, unless noted otherwise. 
In order to investigate the difference between the Wigner 
and streaking delays, we modify the infinite range of the 
potential Eq. (|18[) with a Wood-Saxon factor, 



V a (z) = V c (z) 



1 - 



1 



1 + e~d 2 l~ Zc )/ a 



(21) 



The range of V s (z) is controlled by z c . The parame- 
ter a defines the lengths over which V c {z) is screened to 
approach zero. In our calculation, we use a = 1. For 
z c — > 00 V s converges to V c [Fig. Ufa)]. We numerically 
verified that for values of z c larger than a few atomic 
units, the ground-state wavefunction and energy in V s (z) 
are practically independent of z c [Fig. HJb)]. 

In Fig. [5] and Table |T] we compare the Wigner and 
streaking delays at three different ranges z c . The com- 
parison shows that for an interaction range slightly larger 
than the extent of the ground-state probability distribu- 
tion, say, Zd ~ 2, the Wigner and streaking delays coin- 
cide [see Fig.[5](a)]. In contrast, for z c ^> z^ \tw\ > \ T s\, 
and both delays approach their values for the ID model 
potential Eq. {jXHJ) . In this case the difference of the two 
delays is largest at lower photoelectron kinetic energies 
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FIG. 3: (Color online) Comparison of the Wigner delays r, 
and the streaking delay ts for photoionization of ID model 
hydrogen atoms as a function of the XUV-photon energy fkox ■ 
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FIG. 4: (Color online) The modified Coulomb potential pTf 
(a) and its ground state wavefunction (b) for three different 
interaction ranges z c . 

[Fig. 0(b) and (c)]. Comparison of the three graphs also 
shows that the Wigner delay is more sensitive to changes 
in z c than the streaking delay. 

We emphasize that Wigner delays are calculated with- 
out including the action of an IR-laser pulse, which might 
give rise to the question of whether we should instead 
compare the streaking delay with the Wigner delay ob- 
tained in the same IR field, as mentioned in the intro- 
duction. Obviously, inclusion of the IR-laser electric field 
would change the Wigner delay Q ■ In this case the pho- 
toelectron velocity in Eq. ([7]) would depend on the XUV- 
IR delay At, and the use of this equation would deter- 
mine a Wigner delay that varies with At and the 
IR-pulse intensity. The streaking delay, in contrast, does 
not depend on At. It is also independent on the inten- 
sity of the streaking laser, if this intensity is sufficiently 
low 0, @, [1] ■ Therefore, it is only meaningful to compare 
streaking delays with IR-field-free Wigner delays. 

By computing photocmission spectra and streaking dc- 
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FIG. 5: (Color online) Comparison of the Wigner delay tw 
and the streaking delay ts for three interaction ranges z c in 
Eq. (f2Tjl as a function of the XUV-photon energy huix ■ 
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TABLE I: Streaking and Wigner time delays from Fig. [5] for 
three range parameters z c and three XUV-photon energies 



lays with and without including the IR vector potential 
in Eq. (fl2|). we found that for this atomic target polar- 
ization effects of the initial state of the active electron in 
the electric field of the streaking laser are negligible. This 
is due to the large energy gap between the ground state 
and the excited states ja|. In contrast, the initial state 
polarization is relevant for the case of photoemission from 
solid surfaces discussed in the following section. 



B. One-dimensional Solid 

We model the ID solid surface as a row of N equidis- 
tant atomic layers and represent each atom by a Gaussian 
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FIG. 6: (Color online) (a) Model lattice potential consisting of 
equally-spaced Gaussian potential wells, (b) Corresponding 
band structure with one core-level band and one conduction 
band. The dotted line indicates Fermi level ef- 



potential well to form the lattice potential 

An 



V latt (z) = V - 



N 



2na 



-[z+(i+0.5)a iQtt ]7(2<T ) 



(22) 



where ai a tt is the lattice constant, a controls the overlap 
of the two adjacent atomic potentials, and Vq and A 
are chosen to match the known Fermi energy. We have 
oriented the z-axis with increasing values towards the 
vacuum side and put the origin (z = 0) at the distance 
0.5ai att in front of the top nucleus. Diagonalizing the 
time-independent Schrodinger equation 

2dz 2 



+ Vl a tt(z) 



Mz) (23) 



for N = 47, Vo=-0.5, A Q =2, ai a tt = 6, and a = 0.1o/ att , 
we obtain the core-level and conduction-level bands 
shown in Fig. [6j The Fermi energy is ef = —10.9 eV. 

For the calculation of the XUV photoemission spec- 
trum, we replace V(z) in Eq. ([9]) by Vi a tt{z) and add 
the damping term —iv z /(2X), with the velocity v z = 
y / 2(ujx — \en\), to Viatt(z), in order to model scatter- 
ing of the photoelectrons inside the solid. In a previous 
study we adjusted the electron mean-free path inside the 
solid to A = 5 A as which corresponds to the minimum 
of the universal curve for A as a function of the electron's 
kinetic energy [l4|- We continue to use this value for the 
present investigation. In the calculation of the streaked 
spectra, we further assume an exponential damping of 
the IR-laser field inside the solid, 



A L (z,t) = A L (t) \e z ' 5L Q{-z) + Q{z) 



(24) 
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FIG. 7: (Color online) Streaking and Wigner delays, rs and 
Ty^r , for XUV photoemission with hwx = 100 eV from the 
core-level band of a ID model solid surface as a function of the 
core-level index n. All delays are computed for the electron 
mean-free path A = 5 A and no penetration of the IR-laser 
field into the solid (Sl = 0) for the streaking delay. 

characterized by the IR-skin depth Sl- 

Within each band, we use the index n to label Bloch 
wave functions ip n with energies e„, starting with E\ for 
the lowest energy Bloch wave at the band bottom. Each 
initial Bloch wave below the Fermi level contributes to 
the photoemission spectrum with the energy-differential 
emission probability 



P n = P(e n ,At) 



i|^„(fc,oo;At)| 2 . (25) 



Similarly, each Bloch wave contributes to the COE 
EcoE,n(At), Wigner delays t^^, and streaking delays 
T~S.n- We first calculate the band-averaged COE 

1 



EcoE(At) 



PnEcOE,n{ At ) 



(26) 



£„<£F 



separately for each band. Next we use Eq. (jTTJ) to obtain 
the band-averaged streaking delay rs . Similarly, we find 
the band-averaged Wigner delays according to 



iji 



1 



IJI 
W,n ■ 



(27) 
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A comparison of streaking and the two Wigner delays 
for an XUV-photon energy Hlux = 100 eV and emission 
from the core-level band is shown in Fig. [7|as a function 
of the core-level index n. The monotonic increase in level 
index n in Fig. [7] is just a coincidence. For other ener- 
gies, these delays do not necessary increase with n. As 
for the atomic case (cf., Sec. IIIII A). we find that the dif- 
ference between the two Wigner delays is negligible. We 
therefore only present results for r^J which we denote 
simply as tw below. All numerical results shown below 
are converged in the number of included atomic layers 
N. 

Wigner and streaking delays as a function of the XUV- 
photon energy hwx for photoemission from three indi- 
vidual core- and three conduction-band Bloch levels are 
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FIG. 8: (Color online) Wigner (a) and streaking (b) time de- 
lays for XUV photoemission from three core-level Bloch states 
with energies £5,525, and £47 given relative to the ionization 
limit. Streaking delays are computed for the electron mean- 
free path A = 5 A and no penetration of the IR-laser field 
into the solid (Sl = 0). 
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FIG. 9: (Color online) Wigner (a) and streaking (b) time de- 
lays for XUV photoemission from three conduction-band lev- 
els with energies £i,£io, and £30 given relative to the ioniza- 
tion limit. The occupied part of the conduction band extends 
from £1 to the Fermi level at ef — —8.22 eV. Streaking delays 
are computed for the electron mean-free path A = 5 A and 
no penetration of the IR-laser field into the solid (Sl = 0). 



shown in Figs. [8] and [9l respectively. According to its 
definition Eq. (J7J), the Wigner delay can be regarded as 
an "effective" propagation time for the photoelectron to 
emerge from the solid. For special case Sl = 0, the 
streaking delay is the travel time photoelectrons need be- 
fore getting exposed to the streaking IR field outside the 
solid. Therefore, intuitively, for Sl = only, one would 
expect the streaking delay to be almost identical to the 
Wigner delay. As Figs. [71(101 show, this is confirmed by 
our numerical results. 
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FIG. 10: (Color online) (a) Band-averaged Wigner and 
streaking delays for XUV photoemission from the (a) con- 
duction band and (b) core-level band. Streaking delays are 
computed for the electron mean-free path A = 5 A and no 
penetration of the IR-laser field into the solid [5l = 0). 

It is interesting to observe that, for emission from the 
core-level band, the band-averaged Wigner and streak- 
ing delays decrease monotonously with increasing huix 
[Fig. [10] (a)]. This decrease closely follows the effec- 
tive propagation time \/\/2e of photoelectrons inside 
the solid prior to reaching the solid- vacuum interface at 
z = 0, even though delay contributions from individual 
Bloch levels, r„ (see Fig. [8]), do not show this behav- 
ior. This confirms the interpretation that band-averaged 
Wigner and streaking delays for emission from the core- 
level band can be regarded as an average time needed for 
a released photoelectron to travel a distance A inside the 
solid dp. 

However, this interpretation is not valid for photoemis- 
sion from the conduction band, where the band-averaged 
delays behave non-monotonously as a function of fiiox 
as shown in Fig. fTOf b). A possible explanation for this 
difference is the delocalized nature of the conduction- 
band Bloch wave. This can be checked by examining 
the core-level- band-averaged delay as a function of the 
overlap parameter a in Vi a tt(z) and will be discussed in 
a forthcoming publication [l5[. This means that, even 
for Sl = 0, the relative delays between photoemission 
from the core-level and conduction band can not be due 
solely to the photoelectron's average travel time in the 

solid dini. 

Since the i) equivalence of t$ and tw and ii) the in- 
terpretation of streaking delays in terms of an effective 
propagation time in the solid are only valid for the spe- 
cial case Sl = 0, we next investigate the dependence 
of photoemission delays on the IR-skin depth Sl- Fig- 
ure [TTJ shows the band- averaged streaking delay for emis- 
sion from the core-level band for two XUV-photon ener- 
gies. These results show a very sensitive dependence of 
ts on the IR penetration depth, with t$ changing from 
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FIG. 11: (Color online) Band-averaged streaking delays for 
XUV photoemission with photon energies of hwx = 50 and 
100 eV as a function of the IR-laser skin depth Sl- The elec- 
tron mean-free path is A = 5 A. 
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FIG. 12: (Color online) Band-averaged streaking delay in the 
photoelectron spectrum for photoemission from the conduc- 
tion band for three values of the Fermi energies ef- A = 
5 A and Sl = are used. 

positive to negative delays. In contrast to photoemission 
from the energetically isolated ground states of atoms 
(Sec. MI A[) . the N Bloch waves form a quasi-continuum 
and can be easily hybridized in the IR-laser electric field. 
This initial-state hybridization effect is the stronger the 
deeper the IR electric field penetrates the solid and ac- 
counts for the Sl dependence of rs- We note that the 
actual IR-skin depth is much larger than the electron 
mean-free path, such that the effective depth over which 
photoelectrons are assembled is limited by A, and the 
IR-skin depth tends to become irrelevant for the pho- 



tocurrent Q. Accordingly, the photoemission delays in 
Fig. [TT1 converge in the limit of large #l/A. 

Another interesting observation, as shown in Fig. Q21 is 
that the band-averaged Tw and ts from the conduction 
band depend on the position of the Fermi level. There- 
fore, changing the occupation probability of Bloch waves 
in the conduction band, for example, by increasing the 
temperature or by doping, will change the band-averaged 
delay. This dependence is absent in the core-level band 
because all the core levels are below the Fermi level and 
fully occupied. 

IV. CONCLUSIONS 

We have examined the relation between Wigner and 
streaking delays in the XUV photoelectron emission 
from model atoms and solid surfaces. We showed that 
both, the creation of the photoelectron and its propaga- 
tion contribute to the Wigner delay. For photoemission 
from atoms, the two delays are only identical for short- 
range ionic potentials. For photoemission from surfaces, 
Wigner and streaking delays become identical only in the 
limit of no IR-field penetration into the solid, and both 
delays can be interpreted as the travel time of the photo- 
electron to the surface. Furthermore, for electron emis- 
sion from the core-level band, both delays can be under- 
stood as the average time photoelectrons need to travel 
a distance equal to the mean-free path in the solid. This 
interpretation docs not hold for photoelectron emission 
from the conduction band. This dissimilarity is expected 
to be due to the different (localized versus delocalized) 
nature of core and conduction-band levels. 

We find that streaking delays are very sensitive to 
changes in the IR-skin depth and Fermi energy and de- 
viate from Wigner delays for non-zero IR-skin depths. 
Their dependence on the substrate temperature, impu- 
rities, and adsorbate coverage may leave a measurable 
signature in relative photoemission delays. 
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